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Abstract 

We analyze the neutron electric dipole moment (EDM) in the Minimal Super- 
symmetric Model with explicit i?-parity violating terms. The leading contribution 
to the EDM occurs at the 2-loop level and is dominated by the chromoelectric 
dipole moments of quarks, assuming there is no tree-level mixings between sleptons 
and Higgs bosons or between leptons and gauginos. Based on the experimental 
constraint on the neutron EDM, we set limits on the imaginary parts of complex 
couplings X'ijk and Ajj^ due to the virtual 6-loop or r-loop. 



The minimal supersymmetric standard model (MSSM) has been widely considered 
as a leading candidate for new physics beyond Standard Model. However, unlike the 
standard model, supersymmetry allows renormalizable interactions which break i?-parity 
defined as (—1)3^+-^'+^ and violate the lepton and/or the baryon numbers. It is in fact one 
of the main theoretical weaknesses of these models because -R-parity conservation is an ad 
hoc imposition which may or may not have a fundamental theoretical basis. Therefore, 
instead of neglecting them completely, it is interesting to ask how small these i?-parity 
breaking (/?) couplings could be by investigating directly phenomenological constraints 0|. 

The most general renormalizable /^-violating superpotential using only the MSSM 
superfields is 

Wi^ = K,tUL,E^ + K^uUQjDl + K^^tDpi + • (1) 

Here, are generation indices. The couplings Af^ and Xf- must be antisymmetric 

in flavor, Ajj^ = —\jik and X'lj^. = —X'-^j- There are 36 lepton number non- conserving 
couplings (9 of the A type and 27 of the A' type) and 9 baryon number non- conserving 
couplings (all of the A" type) in Eq.(l). To avoid rapid proton decay, it is usually assumed 
in the literature that A, A' type couplings do not coexist with A" type couplings. This 
can be achieved easily by imposing baryon number symmetry. The bilinear terms fijLjH2 
contribute to lepton flavor and number violation and could be responsible for neutrino 
masses. Phenomenologically, many of these couplings have been severely constrained 
using low-energy processes or using high energy data at the colliders 0-|jTO[- In this 
paper, we shall not consider A^'^^ and fij couplings. 

However, most of the bounds in the literature constrain the real part of the trilinear 
couplings, or the product of trilinear couplings. The exception is the bound coming from 
the ex which constrains Im(A^^2-^i2i) < 8 ^ 10^^^ [O]. We propose to study the neutron 



electric dipole moment, which is tightly bound by experiment, and thus obtain limits on 
the imaginary parts of different products of trilinear couplings from the ones imposed by 
ex- 

The electric dipole moment of an elementary fermion is deflned through its electro- 
magnetic form factor F^{q^) in the (current) matrix element: 

{f{p')\JM\f{p)) = ^WM)<p)^ (2) 
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where q = p' — p and 

r^l?) = i^i(g')7/. + ^^I^^^M-^' + ^A(g')(7M75g' - 2m^,q^) + ^^a^^^.q'', (3) 

with m the mass of the fermion and -Fi(O) = Cf. The electric dipole moment (EDM) of 
the fermion field / is then given by 

<i/ = ^. (4) 

corresponding to the effective dipole interaction 

Cedm = -'-djfa.^j.fF'^'' . (5) 

In the static limit this corresponds to an effective Lagrangian £edm —>■ df'^^cx ■ E'^a, 
where is the large component of the Dirac field. Similarly the quark chromoelectric 
dipole moment (CEDM) is the coefficient in the effective operator: 



CcEBM = -'-d',qa,u^5Y9G'""' . (6) 



The relevant Lagrangian for generating an EDM is: 
^ = -{lll ^^^^^^^ + H-c. j + ■ ■ ■ = [K,,{-i^u,dl + i)d,dl) + H.C.] + ■ ■ ■ . (7) 

It has been shown that there is no one-loop contributions to EDMs based on A, 
A' or A" couplings based on helicity properties and symmetry. Here we briefly review its 
origin. It is easy to show that one cannot induce EDMs from the diagram that requires 
the external mass insertion due to the equation of motion. As a result of this lemma, 
proper helicities for external fermion lines have to come directly from vertices. Let us look 
at the electron EDM, which needs external L and E^. For the correct quantum number, 
possible one loop contributions have to be proportional to either (1)AA* or (2)A'A'*. Based 
on the above lemma, the external L and E'^ are required to come directly from vertices. 
Case (1) cannot produce the helicity flip. Case (2) is even worse, there is no vertex to 
give -E^. So the one-loop electron EDM is absent. For the d quark EDM, possibilities 
are either (1) A'A'* or (2) A"A"*. Case (1) does not work because both di and d'^ have 
to come from a CP-even product of a complex conjugated pair of vertices, and case (2) 
fails badly because there is no vertex to give an external di. Similar reason follows for 
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the u quark EDM. As a reminder, there are one-loop EDM amphtudes |]12[ related to the 
bilinear term HjLjH2, which mixes sleptons and Higgs bosons etc. We do not consider 
these couplings /ij in this work. 

At the two-loop level, a number of different types of configurations contribute, which 
we classify as rainbow-like (I), overlapping (II), tent-like (III) and Barr-Zee (IV)||1^ types. 
The rainbow-like graphs (I) are those with two concentric boson loops, the outer of which 
must be a charged Higgs loop (for the same reason that 1-loop graphs do not exist). The 
inner loop may be a left or right sfermion. The complete set of this type of graphs is 
given in Fig. la. The complete set of overlapping type of graphs is given in Fig. lb. In 
this case, one of them must be a charged Higgs, the other a left or right sfermion. 
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Fig. la (i) Rainbow-like diagram for the d quark. The generic /R vertex is marked by o 
and its complex conjugate by •. 
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Fig. la (ii) Rainbow-like diagram for the u quark. 
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Fig. lb (i) Overlapping diagram for the d quark using A'. 
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Fig. lb (ii) Overlapping diagram for the u quark using A'. 



5 



* k.R 



X H- 



\ 



dj, 



kR 



UjL S dR 

7 



f 

I 

' 



I HL \ 



dr d, 



kR 



^iL dn 

7 



Fig. Ic (i) Tent-like diagram for the d quark using A'. 
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Fig. Ic (ii) Tent-like diagram for the u quark using A'. 
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Fig. Id Barr-Zee type graph for u quark EDM 
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The tent-like graphs (III) have a trihnear bosonic vertex. Again, the 3 different boson 
legs can be two sfermions and one charged Higgs (in all possible configurations). The 
complete set of this type of graphs is given in Fig. Ic. Careful consideration of all the 
type I-III graphs shows that their contributions are suppressed by both one power of light 
quark mass plus some CKM mixing angles factor compared to those of type IV (Fig. 2). 
Therefore, one expects the Barr-Zee type of contributions dominate and we shall study 
them in detail next. 




d-R diiq) di 



Fig. 2 A typical two-loop diagram of the Barr-Zee type. Note that there are 3 ways to 
insert mass. 

In Ref. 0, only a rough estimate of the two-loop contributions to EDMs is provided. 
We shall present here a complete calculation of the quark (or electron) EDM and the 
quark CEDM at the two-loop level due to the Barr-Zee type mechanism and show that 
the neutron EDM is dominated by the CEDM of the d quark. This calculation leads to 
more stringent bounds than previously obtained. 

There is another class of Barr-Zee graphs with sneutrino line replaced by the charged 
slepton line and corresponding modifications of the fermions charges in the loop. The 
calculations of these type of graphs are very similar to the one in the charged Higgs 



models of CP violation as in Ref.O, 15]. Comparing the charged Higgs contributions 



to the EDM in Ref. p^, [T5| with the neutral Higgs contributions given in Ref. ||13||, one 
can observe that the neutral Higgs contributions generally dominates given comparable 
coupling constants and boson masses. Therefore, we shall only give details of sneutrino 
contributions here. 
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i>gg vertex of the inner loop 



The two-loop diagram to the CEDM of the d-type quark can appear with the coupUng 
X'^jf. through the virtual vertex ugg. The amplitude of the inner loop in terms of the 
leading gauge invariant terms is: 

T'^'' = S(q') [k'^q'^ - k ■ qg^ + P{q^) [ie^^^^Paqp] , (8) 
where S and P correspond to scalar and pseudoscalar form factors respectively: 



IGtt^ Jo ml — y{l — y)q^ ' 



IGtt^ Jo ml — y{l — y)q'^ 

Second loop 

Combining the two twisted diagrams and the two choices of sneutrino flow directions, we 
have a combinatoric factor of 4 in the the two loop CEDM amplitude. In the convention 

of Eqs.(2-6), we obtain the CEDM of the d quark at the scale of mf,, 



Js)^^~ 1287r^ '"^^^^^^^^^^Vo '^io {ml + y{l-y)Q-){Q^ + Ml)Q^ 



327r3 Ml \ Ml 

with the loop function 



+ 2 + lnr+ (lnT)2 for r ^ . (11) 



3 

Implicit sum over sneutrino flavors i is assumed in the above. The last asymptotic form 
is useful because the ratio r = m^/Mj^ is small. The large logarithmic factor helps place 
a strong constraint on A' couplings. Note that sneutrino is the heaviest particle in the 
loop. At mb scale, the sneutrino induces a four fermion interaction of b and d quarks. As 
a result, by simple power counting, the gluonic loop is logarithmically divergent which 
explains the large logarithmic enhancement factor. 
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Replacing the gluon line by the photon line, we obtain the EDM of the quark simply 
by substituting the color factor and the charge factor. 



a. 



(12) 



Now we address the QCD evolution of these Wilson's coefficients. As an effective 
theory, the 4-fermion vertices of the form (hh) (dd) arise first in the energy scale below M^. 
when z> are integrated out of the theory. Then the EDM and the CEDM of the d quark 
arise below scale. Therefore, the A' couphngs in the above equations are evaluated 
at the rrib scale. We ignore the dressing of these 4-fermion vertices because of the small 
value of their couplings and the slow running of Og at such high energy scale. In this 
perspective, the A' factors in Eq. (^ are defined at the short distance scale near Mo. 
Below rrih, the CEDM and the EDM of light quarks appear and they evolve down to the 
hadronic scale Ah by 
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(13) 



(14) 



Note that in some references |1^, a light quark mass coefficient has been factored out so 
that the form of evolution equation looks different from above. We denote by D!^ {DX) 



the neutron EDM due to the CEDM (EDM) of light quarks. The SU{Q) relation gives: 

(15) 
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For a,{Mz) = 0.12 and 5(,(AH)/(47r) = 1/V6, the QCD evolution factors Z'^ and Z^ are 
about 0.71 and 0.84 respectively. Our formulas and numerical values are consistent with 
those in Ref. but differ from those in Ref. JTB|. 

For completeness, we add another large contribution to the d quark EDM due to 
the r-lepton replacing the b quark inside the first loop. We obtain two independent 
contributions as: 



^ / mi, i=l,2,3 



Ml 



Im(AS3A:n) ■ F 



(16) 
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The latter contribution from the r-loop is induced at the scale and we need to adjust 
the minor change in the QCD evolution. There are also other Barr-Zee type diagrams 
from the exchange of or Z gauge bosons. However, they are known to be giving 



smaller contributions and thus we ignore them in our numerical study |T5|, |I6|. 

As is not directly involved in the fi interaction, the u quark CEDM do not appear 
through A' in the form of Fig. 2. Nonetheless, there are two- loop diagrams Fig. la (ii), 
Ib(ii), lc(ii) and Id which are suppressed by the light quark mass and mixing angles. 
Therefore, the fi contribution to the neutron EDM is dominated by the d qaurk CEDM 
and EDM. Assuming all My. are equal and taking typical values, M^^ ^ 300 GeV and 
rrib ~ 4.5 GeV, we have 

Di ^ 5.46 X 10-21 (e-cm) x MKssKn) , (18) 

i 

Dl ^ -1.03 X 10-22 (e-cm) x ^ \ra{K^,K^,) 

i 

-1.92 X 10-22 (e-cm) x ^ Im(A*33A:n) ■ (19) 

Our numerical result shows that the strongest constraint comes from the CEDM of the 
d quark. Using the up-to-dated experimental [jl9| bound < 6.3 x 10-26 e-cm and 
barring accidental cancellation among contributions, we derive the constraints: 

< 1.2xl0-^ (20) 

i 

i:im(A*33A:,,) < 33X10-S (21) 

for My = 300 GeV. 

In Fig. 3 we plot both the photon and gluon contributions to the neutron EDM versus 
the sneutrino mass My in the region of interest (100 to 600 GeV) with J^i^^T^i^'issKii) 
or X^i Ini('^i33'^iii) scaled to 10-^. One could see that the gluon contribution consistently 
dominates the corresponding photon one by at least an order of magnitude over the whole 
parameter space explored. 
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Fig. 3 The neutron EDM Dn versus Mjj with J2i'^^i>^'i33Kn) or J^i'^MKssKii) scaled to 10"^ 



The electron EDM can arise via both A' or A-type -R-parity violating coupling Ref . |^ . 
Based on above study, the analytical formula for the electron EDM at the two-loop level 
is: 



E 



Ml 



ill. 



■F 



+ ]^Im(A;33A2n) • F 



[mi 

W 



[22] 



In Fig. 4 we assume all Mj>^ to be equal and plot contributions to the electron EDM versus 
the sneutrino mass My in the region of interest (100 to 600 GeV). Using the up-to-dated 
experimental |2y] bound |(ie| < 0.43 x 10~^^ e-cm and barring from accidental cancellation 
among contributions, we derive constraints: 

Im(A;33A2ii) < 0.74 x 10"^ , (23) 

$:im(Ar33Aai) < 1.3xl0-\ (24) 
i^i 



for Mo = 300 GeV. 
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Fig. 4 The electron EDM de versus Mp. 

Conclusion 

We have presented an exact and complete calculation of the dominant contribution to the 
neutron EDM in a minimal supersymmetric model without R parity due to the couplings 
A and A'. The CP violation does not depend on the complex phases 0^ and (pAo (the 
phases of the Higgsino mass parameter and the trilincar scalar coupling Aq) in Minimal 
Supergravity models, therefore unrelated to the restrictive bounds or complicated can- 
cellations necessary in MSSM. The leading contribution to the neutron EDM occurs 
at two-loop level through the Barr-Zee mechanism. We obtain stringent bounds on the 
product ImAggA'iii < 0(10-^). 
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Figure Captions 

Fig. 1 a. (i) Rainbow-like diagram for the d quark. The generic vertex is marked by 
o and its complex conjugate by •. (ii) Rainbow-like diagram for the u quark. 

b. (i) Overlapping diagram for the d quark using A', (ii) Overlapping diagram for 
the u quark using A'. 

c. (i) Tent-like diagram for the d quark using A', (ii) Tent-like diagram for the u 
quark using A'. 

d. Barr-Zee type graph for u quark EDM 

Fig. 2 A typical two-loop diagram of the Barr-Zee type. Note that there are 3 ways to 
insert mass. 

Fig. 3 The neutron EDM D^^ versus Mp with ^M^'^sKu) or Ei Im(A*33A^ii) scaled to 
10-^ 

Fig. 4 The electron EDM de versus Mj>. 
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